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Abstract 



We discuss the parsec-scale structural variability of the extragalactic jet 3C 111 related to 
a major radio flux density outburst in 2007. The data analyzed were taken within the scope 
of the MOJAVE, UMRAO, and F-GAMMA programs, which monitor a large sample of the 
^ . radio brightest compact extragalactic jets with the VLBA, the University of Michigan 26 m, 

Y'{ ' the Effelsberg 100 m, and the IRAM 30 m radio telescopes. The analysis of the VLBA data is 

C/3 , performed by fitting Gaussian model components in the visibility domain. We associate the 

ejection of bright features in the radio jet with a major flux-density outburst in 2007. The 
evolution of these features suggests the formation of a leading component and multiple trailing 
, components. 
^ ' Keywords: galaxies: individual: 3C 111 - galaxies: active - galaxies: jets - galaxies: nuclei 

a^ . 

1 Introduction 

Jets of active galactic nuclei (AGN) are among the most fascinating objects in the universe. From 



the time when the term "jet" was first introduced by Baade & Minkowski (1954) until today it 
is still unclear how these jets are created and formed. A prime source to gain insight into the 
physics of extragalactic jets is the broad-line radio galaxy 3C 111 (PKS B0415-I-379) at z = 0.04!|3. 
The object can be described with a classical FR II morphology ( [Fanaroff &: Riley 1974[ ) exhibit- 
\rl \ ing two radio lobes with hot spots and a single-sided jet ( [Linfield fc Perley 1984[ ). Untypical for 

. radio galaxies, a small inclination angle of only 18 ° to our line of sight has been determined on 

\ parsec scales ([Jorstad et al. 2005| . Moreover, 3C 111 has a blazar-like spectral energy distribution 

(SED) (jHartman et al. 2008^ and shows one of the brightest radio cores in the mm-cm wave- 
length regime of all FRII radio galaxies. Superluminal motion was detected in this radio galaxy 



by Gotz et al. (1987) and Preuss et al. (1988) making this source one of the first radio galaxies 
to exhibit this effect. The EGRET source 3EG J0416-I-3650 has been associated with 3C111 
([Sguera et al. 2005 IHartman et al. 2008P and 7-ray emission from 3C 111 has been confirmed by 



Fermi/LAT (|Abdo et al. 2010a[Ebdo et al. 2010bllAbdo et al. 2010cp . A major flux-density out- 



burst in 1997 was investigated by Kadler et al. (2008) with 10 years of radio monitoring data 



(1995-2005). In addition Chatterjee et al. (2011) and Tombesi et al. (2011) report on a possible 



connection between the accretion disk and the jet of 3C 111. 

In this paper a new major flux density outburst from 2007 and the associated jet kinematics will 
be discussed with data from the Very Large Baseline Array ( VLBA0) , the University of Michigan 
Radio Astronomy Observatory (UMRACB) and the F-GAMMA program. 



1 Assuming Ho = Tlkms-^ Mpc-^.^A = 0.73 and = 0.27 (1 mas = 0.95 pc; Imasyr^^ = 3.1c) 

•^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 

cooperative agreement by Associated Universities, Inc. 

''UMRAO has been supported by a series of grants from the NSF and NASA and by funds from the University 

of Michigan 
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Figure 1: Long-term radio lightcurve of 3C 111 obtained by the UMRAO at 4.8 GHz (left). Short- 
term radio hghtcurves (right) at 32.0 GHz, 42.0 GHz, 86.24 GHz and 142.33 GHz obtained by the 
F-GAMMA program. 



2 Data Analysis 

The broad-hne radio galaxy 3C 111 has been part of the VLBA 2cm Survey program since its start 
in 1995 (jKellermann et al. 1998[) and its successor MO JAVE (Monitoring Of Jets in Active galactic 
nuclei with VLBA Experiments, ILister et al. 2009| in 2002. Twenty-four epochs of data have been 
taken from 2006 to 2010 within the MOJAVE program of this source. Phase and amplitude self- 
calibration as well as hybrid mapping by deconvolution techniques were performed as described 



by Kellermann et al. (1998) Utilizing the program DIFMAP ( [Shepherd 1997[ ), two-dimensional 
Gaussian components have been fitted in the (M,'!;)-plane to the fully calibrated data of each epoch. 
We refer to the inner '-^ 0.5 mas as the "core" region, which can usually be modeled with two Gaus- 
sian components. All models have been aligned by assuming the westernmost components to be sta- 
tionary and all component positions are measured with respect to it. Conservative errors of 15% are 
assumed for the flux densities of the model-fit components accounting for absolute calibration un- 
certainties and formal model- fitting uncertainties (jHoman et al. 2002]) . Within the UMRAO radio- 
flux-density monitoring program (jAller et al. 2003p . more than three decades of single-dish flux- 
density data have been collected for 3C111 at 4.8 GHz, 8.0 GHz, and 14.5 GH^ In addition, 3C 111 
is observed monthly by the F-GAMMA program (jFuhrmann et al. 20071 [Angelakis et al. 2008| at 
multiple frequencies throughout the cm- and mm-bands since 2007. 



3 Results 

3.1 Lightcurves 

Figure [1] (left) shows the long-term radio lightcurve of 3C 111 at 4.8 GHz. Since the start of the 
measurements the source has been in a low activity state for almost two decades with only minor 
activity. The source showed a major outburst in 1996/1997 (!K adler et al. 2008|) . Starting in 2004, 
minor outbursts are observed and the flux density level increases. A major flux density outburst 
starts in early 2007 at high frequencies (see Figure [U right), peaking ~ 2007.6 and is subsequently 
seen at lower frequencies. A secondary outburst starts mid 2008 at high frequencies. The overall 
flux density level at 4.8 GHz has been decreasing since the major outburst. 

3.2 VLBA Data Overview 

In this work, we focus on the time period 2006 through 2010, which contains 24 MOJAVE epochs 
with an average image rms noise of 0.2mJybeam~^ and a maximum of 0.4mJy beam~^. The 
restoring beams of different epochs are very similar with an average of (0.84 x 0.58) mas at P.A. 
—8°. A peak flux density of 6.54 Jy was measured in May 2008. Figure [2] shows the evolution 
of the parsec scale jet of 3G111 observed within the MOJAVE VLBA program at 15 GHz since 

*In this work, we consider only the 4.8 GHz UMRAO data. A multi-frequency long-term analysis of the light 
curve will be presented elsewhere (Grossberger et al., in prep.). 



2 




2010-12-24 



Figure 2: Naturally weighted clean images of 3C 111 from 2008 to 2010. The minimum contours 
are set to 5 sigma of the rms noise and increase logarithmically by a factor of 2. All fitted model 
components are indicated with a circle or ellipse (size of the full-width half maximum of the 
Gaussian function) enclosing a cross. The lines are the fitted evolutional tracks of the leading 
components. The dashed line indicates the primary component (and leading component after mid 
2009). The dotted line shows the position of the core. 
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May 2008 (excluding six epochs before the ejection of the primary component). A major bright 
feature appears in May 2008 and is travelhng downstream. The source brightness distribution was 
modelled with the CLEAN algorithm by|Hogbom (1974)| within DIFMAP. 



3.3 Model Fitting 

In Figure [31 the radial distances to the stationary component in the core region are plotted for 
every component as a function of time. The identification throughout the epochs is based on the 
comparison of the positions and flux densities of these model components. This identification is 
preliminary and will be discussed in more detail in a forthcoming paper. In the context of this 
paper we focus on the components which can be associated with the major flux density outburst 
of 2007 (see Fig. H]). 

A linear regression fit has been performed to measure component speeds and ejection dates (see 
Figure [3|). The ejection dates of a primary component and a secondary component are quasi iden- 
tical (~ 2007.6) within the errors and were found to coincide with the peak of the outburst at high 
frequencies (see Figure [U right). The components flux density evolution (see Figure |4]) shows that 
the core region was extremely bright during the time of the outburst but dropped significantly 
after ejecting the bright primary and secondary component. The determined apparent speed for 
the primary component is 3.94 ± 0.19 c and for the secondary component 2.80 ± 0.40 c. 
The primary component remains at a constant flux density of IJy until mid 2009. After that, 
the component is splitting into multiple parts: a new leading component with trailing components 
in its wake. 

The secondary component has a higher flux density than the primary in the beginning of its lifetime 
which rapidly decays. This decay suggests that this component disappeared in mid 2009 though 
an identification with the first trailing component is possible based on position alone. The calcu- 
lated apparent speed of the new leading component after mid 2009 is 4.53 ± 0.09 c with the flux 
density decreasing. The first trailing component has an apparent speed of 3.32 ± 0.19 c and shows 
a constant flux density evolution. The second trailing component was flrst observed 2010.53 with 
a flux density of ~ 240 mJy and could be modeled until 2010.91 with a flux density of less than 
lOOmJy. The flux density evolution of the second trailing component suggests that this component 
faded away and thus could not be modeled in epoch 2010.98. An association of this second trailing 
component with the 2008 component is possible based on the flux density evolution and position 
but needs further investigation of the 2008-outburst. 

A similiar behaviour with a leading, secondary and trailing components has been seen in the evo- 



lution of the components associated with the outburst from 1997 by Kadler et al. (2008) The 
components were interpreted as a forward and backward shock with the backward shock fading 
very fast. In this model, the plasma of the forward shock entered a region of rapidly decreasing ex- 
ternal pressure allowing it to expand into the jet ambient medium and accelerate. In the following, 
the plasma recoUimated and trailing features were formed in the wake of the leading component 
(jPerucho et al. 2008p . 



4 Summary 

In this paper, the ejection of new jet components on parsec scales were associated with a major 
flux density outburst of 3C 111 in 2007. It was shown that the major flux density outburst can be 
associated with the ejection of a primary jet component and secondary component. The evolution of 
the leading component suggests a split into multiple components. The full multi-epoch kinematical 
analysis of the VLBA jet of 3C 111 between 2006 and 2011 will be presented elsewhere (Grossberger 
et al., in prep). 

Abdo, A. A., et al., 2010a, ApJ 715, 429 
Abdo, A. A., et al., 2010b, ApJS 188, 405 
Abdo, A. A., et al., 2010c, ApJ 720, 912 

AUer, M. F., AUer, H. D. & Hughes, P. A., 2003, in Radio Astronomy at the Fringe, Zensus, J. A., 

Cohen, M. H., Ros, E. (eds.), A SP Conf. Ser. 300, 159 
Angelakis, E., et al., 2008. larXiv:0809.3912i 



4 




u N -I I I I I I I I -r~ii-~Jr^i u I- I I iX 

2006 2007 2008 2009 2010 2011 2007 2008 2009 2010 2011 



Epoch |y] Epoch [y] 

Figure 3: Separation of all model-fit components (left) and only components associated with the 
2007-outburst (right) with respect to the core as a function of time. Model-fit components which 
could not be identified in 5 or more epochs are marked with a black cross. The 2007-primary com- 
ponent is presented by triangles, the 2007-secondary component by crosses, the leading component 
by open rectangles, the first trailing component by open diamonds, the second trailing component 
by pluses and the 2008-component by stars. Linear regression fits determine the trajectories of the 
components associated with the 2007-outburst. 
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Figure 4: Flux density evolution of the "core" region and the components as a function of time. 
Conservative errors of 15% are assumed for the flux densities of the model-fit components. The 
component symbols are the same as in Fig. [3] with the addition of the "core" region, presented by 
stars. The black arrow at the bottom indicates the calculated and almost identical ejection date 
of the primary and secondary component based on the derived jet kinematics. 
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